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We report on a search for the rare decay B° — > -D*°7, which in the standard model is dominated 
by W-exchange. The analysis is based on a data sample comprising 87.8 x 10 6 BB pairs collected 
with the BABAR detector at the PEP-II asymmetric-energy B Factory at SLAC. No significant 
signal is observed, and an upper limit on the branching fraction of 2.5 x 10~ 5 at the 90% confidence 
level is obtained. 



PACS numbers: 12.39.St, 13. 20. He 

Within the standard model (SM) , the rare decay B° — > 
£>*°7 is dominated by the VF-boson exchange pro- 
cess. One of the leading SM contributions to the de- 
cay is illustrated in Fig. ^ Similar W-cxchangc tran- 
sitions are present in other decays. For example, they 
contribute to the decay B° — > p°j along with the lead- 
ing electromagnetic-penguin process |2j. The branch- 
ing fraction B(B° — > D*°j) is estimated to be of order 
10~ 6 0, , but the presence of a large qqg (color octet) 
component in the wave function of the B meson may re- 
duce the color-suppression enough to raise the branch- 
ing fraction by a factor of about 10 Q. A search for 
B° -»■ D*°7, published by the CLEO collaboration 0, 
resulted in a limit of B(B° -> D*°~f) < 5.0 x 10" 5 at the 
90% confidence level (C.L.). 




FIG. 1: W-exchange is the leading contribution to the B° — ^ 
D*°7 decay in the standard model. The photon may be emit- 
ted from any quark line or the W. 

We search for the decay B° — ► D*°7 in data collected 
using the BABAR detector operating at the Stanford 
Linear Accelerator Center (SLAC) PEP-II asymmetric- 
energy e + e~ collider. The collider runs with a center-of- 
mass (CM) energy of 10.58 GcV at the peak of the r(4S) 
resonance, which decays into B + B~ and B°B° pairs. The 
analysis is based on 87.8 x 10 6 BB pairs, corresponding 
to an integrated luminosity of 79.9 fb _1 . The BABAR de- 
tector is described in detail in Ref . ; here we introduce 
briefly the detector systems important for the present 
analysis. Tracks of charged particles and their momenta 
are measured in a vertex tracker, consisting of five layers 
of double-sided silicon microstrip detectors, and a 40- 



layer drift chamber. Both systems are located within a 
1.5-T solenoidal magnetic field and provide dE/dx mea- 
surements for particle identification (PID). A Chercnkov 
ring imaging detector adds measurements for PID by 
recording Chercnkov light emitted from charged parti- 
cles traversing transparent quartz bars. Photons are 
identified by an electromagnetic calorimeter consisting 
of 6580 CsI(Tl) crystals. 

Event samples from Monte Carlo (MC) simulations are 
used to optimize the event selection criteria and to esti- 
mate the signal efficiency and background. The detector 
response is simulated using GEANT4 Q. The MC sam- 
ple for the signal B° — > B*°7 contains 328 000 events. 
We use MC samples of similar size for several exclu- 
sive B-decay background modes. The color-suppressed 
hadronic decay B° — > D*°tt°, with branching fraction 
(2.7 ± 0.5) x 10~ 4 [3, is the largest contributor among 
them. Other backgrounds originate from BB modes with 
incompletely or incorrectly reconstructed particles, and 
from random combinations of particles from two different 
B mesons or from qq pairs. For these, we use MC samples 
of generic BB events and continuum qq (q = u, d, s, c) 
events corresponding to about 200 fb -1 and llOfb -1 , re- 
spectively. 

The D*° candidates are reconstructed in six sub- 
modes, with D*° -> B (7T°,7) and D° -> 
(K~ir + , K~tt + it , A'~7r + 7r + 7r~). The event selection cri- 
teria are optimized by using the MC samples to maximize 
S 2 /(S + B), where S (B) is the number of signal (back- 
ground) events. A signal branching fraction of 10~ 6 is 
assumed during the optimization. The most important 
selection requirements are described below. 

The photon from the decay B° — > B*°7 is emit- 
ted with an energy of about 2.3 GeV in the CM frame 
("hard photon"). Although this high energy leads to 
a relatively clear signal, care must be taken that rem- 
nants of 7T° decays are not mistaken as the signal pho- 
ton. The "7T° veto" rejects a hard photon candidate if 
its combination with any other photon with laboratory 
energy larger than 30 MeV yields an invariant mass in 
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the range [110, 155] McV/c 2 . A similar veto for r\ decays 
rejects a photon candidate if its combination with any 
other photon of laboratory energy larger than 250 MeV 
yields an invariant mass within [508, 588] MeV/c 2 . Hard 
photon candidates must also pass a calorimeter shower- 
shape requirement designed to exclude irregularly shaped 
showers caused, for example, by overlapping photons 
from 7T° decay. Background is further suppressed by 
requiring a hard photon candidate to be isolated from 
all other showers and tracks by at least 50 cm in the 
calorimeter. 

A photon candidate from the decay D*° — ^ D°j ("soft 
photon" ) must satisfy the same shower-shape require- 
ment and rj veto that are applied to hard photons. In the 
7T° veto the minimum energy for the other photon is 
raised to 80 MeV and the invariant mass range is re- 
stricted to [115, 150] MeV/c 2 . In addition, the CM energy 
of the soft photon candidate has to be at least 110 MeV. 

The mass of the n° in the decay D*° — > £>°7r° and 
of the 7T° in the decay D° — ► K~ir + ir° is required to be 
within 11 McV/c 2 of the true ir° mass (which corresponds 
to a cut at about 1.7 a, where a is the 7r° mass resolu- 
tion). Photons from ir° decay need a minimum energy 
of 30 MeV and have to pass a similar, but slightly less 
stringent, shower-shape requirement as the hard and soft 
photons. 

The charged K and n tracks are required to originate 
from the interaction point and have to pass likelihood- 
based particle identification selections using dE/dx and 
Cherenkov light measurements. The K track in the 
D — > K~ir + ir + ir~ decay is in addition required to 
have a transverse momentum larger than 0.1 GeV/c and 
at least 12 hits in the drift chamber. A vertex fit 
is applied to the D° candidates. They are required 
to have masses close to the known D° mass: within 
12MeV/c 2 (~ 1.8 a) for D° -> K~tt + , within 23MeV/c 2 
(~ 1.9 a) for D° -> A^tt+tt , and within 12MeV/c 2 
(~ 2.3 er) for D° — > K~ir + ir + ir~ . Additional selection 
requirements are applied to D° candidates decaying into 
K~ir + ir°. The laboratory energy of the ir° must be at 
least 250 MeV, and only D° —> A^tt+tt candidates that 
appear in the Dalitz plot close to known resonances 
are accepted. The difference between the D*° and D° 
mass has to be within 2 MeV/c 2 (~ 2 er) for D*° -> D°tt° 
and within 9 MeV/c 2 (~ 1.8 a) for D*° -► D°j of the 
known value of Ref. || . 

The D*° hclicity angle 0* H is defined in the D* a CM 
frame as the angle between the direction of the D° and 
the direction opposite to the B momentum. For the 
D*° — ► D°n° modes, cos6>^ is distributed as sin 2 0* H for 
signal, but as cos 2 9* H for background from B° — > D*°ir°. 
Optimization leads to the requirement | cos | < 0.75. 
No such condition is imposed for D*° — ► D°j modes. 

Several selection requirements reduce the number of 
fake decays from qq continuum background. The angle 
9* B is defined as the angle between the B candidate mo- 



mentum in the T(AS) CM frame and the beam axis. In qq 
background events the distribution is uniform in cos 9* B , 
while for real B mesons it follows a sin 2 0* B distribution. 
We require that | cos0* B \ < 0.8. The angle 0% is the angle 
between the thrust direction of the B candidate and the 
thrust direction computed from the other photons and 
tracks in the event. For signal events the distribution of 
|cos#y| is flat, while for continuum events the distribu- 
tion has a maximum at | cos#^| = 1 due to their jetlike 
nature. We require that | cos6^| < 0.75. 

The candidates are subsequently characterized with 
two kinematic quantities, mss and AE. For the "energy- 
substituted mass" mES; the energy of the B candidate is 
substituted by precisely known beam parameters: 

™es = y/(s/2 + C 2 Po - PiJ ) 2 /£o 2 -c 2 P 2 B , (1) 

where s is the square of the total CM energy, Eq and po 
are the energy and momentum of the initial T(AS) in the 
laboratory frame, and ps = Pu*o +p 7 is the momentum 
of the B candidate, also taken in the laboratory frame. 
The quantity AE is defined as the difference between the 
energy of the B candidate E* and the beam energy, both 
taken in the CM system: 

AE = E*- l -^7s. (2) 

Requirements of |AA| < 0.34 GeV and 5.2 < m E s < 
5.29 GeV/c 2 are applied at this point. 

If an event contains more than one B° — » D*°j candi- 
date passing all selection criteria, the selection is made 
based on a x 2 function that uses the measured D° mass 
and D*°-D° mass difference, the measured resolutions, 
and known mass and mass-difference values from Ref. . 
This selection is sufficient, as the ambiguity is never due 
to the presence of two hard photon candidates. 

The distribution of tubs versus AE is shown in Fig. |3 
for the data taken at the T(4S) resonance. While 
the combinatorial qq background is smoothly distributed 
over this plane, the signal should peak around AE = 
and rn.Es = 5.28GeV/c 2 . The borders of the signal box 
are given by 5.275 < m E s < 5.285 GeV/c 2 and -0.1 < 
AE < 0.08 GeV, extending to about 1.7 (1.9) times the 
resolution of tues (AE) of signal events. The AE con- 
straint is asymmetric to account for the energy leakage 
from the calorimeter for the hard photon candidates. The 
area with ities ranging from 5.2 GeV/c 2 to 5.27 GeV/c 2 is 
called the "grand sideband." 

The contributions to the systematic uncertainties in 
the signal reconstruction efficiencies are listed in Ta- 
ble^ The overall relative uncertainties range from 16.5% 
to 19.8%, depending on the reconstruction mode (see Ta- 
ble HTJ. The major contributors are described here in 
more detail. The uncertainties in the photon reconstruc- 
tion due to efficiency, energy scale, and energy resolution 
uncertainties are studied with control samples and result 




TABLE I: Maximal and minimal 
tainties in the efficiency for the 
modes. 



relative systematic uncer- 
individual reconstruction 



FIG. 2: Distribution of data events in the AE-itles plane. 
The lines indicate the regions of the signal box and of the 
grand sideband. 



in an uncertainty of 2.5% per photon (5% per tt°). Stud- 
ies of the track finding efficiency using control samples 
result in uncertainties of 2.6% to 5.9% depending on the 
mode. The size of the uncertainty in the AE and mES 
selection is obtained by varying the selection according 
to observed differences between data and MC simulation. 
For the thrust angle the B° angle 6* B , and the heli- 
city angle 6* H , the size of the uncertainties is obtained by 
shifting the selection requirement by ±0.05 in the cosine 
of each angle. The uncertainty due to possible discrep- 
ancies between data and MC simulation in the D° mass 
and the D*°-D° mass difference is estimated by compar- 
ing these distributions for events in the grand sideband. 
Data and Monte Carlo simulation agree sufficiently well, 
and the size of the systematic uncertainty in the efficiency 
is obtained from the uncertainty on the fits to the mass 
and mass-difference plots. 

Several correction factors are applied to the signal ef- 
ficiency based on comparison studies on data and Monte 
Carlo simulations: a tracking efficiency factor of 0.992 for 
the kaon in the decay D° — > K~ir + TT + Tr~ , a factor 0.95 
for the decay D° — > A'~7r + 7r° due to the selection require- 
ment involving the Dalitz structure, and factors from 0.89 
to 0.95 depending on the reconstructed mode due to pho- 
ton reconstruction. The overall selection efficiencies for 
the six signal modes are listed in Table EH The uncer- 
tainties on the efficiencies include all contributions from 
systematic effects on the efficiencies. The combined effi- 
ciency (weighted by the branching fractions of the indi- 
vidual modes and taking correlations in the uncertainties 
between the six submodes into account) is (1.8 ± 0.3)%. 
In the determination of the B° — > D*°j branching frac- 
tion results, a 1.1% uncertainty on the number of BB 
pairs in the data sample is included as well as the con- 
tribution by the D° (D* a ) branching fraction uncertain- 
ties @. 

The number of events expected in the signal box due 
to background is not estimated from data, but from MC 



Systematic uncertainty 
in % of the efficiency 


7 & 7T° reconstruction 


5.0 to 12.5 


Hard 7 separation 


2.0 


Shower shape 


1.0 to 2.5 


7T°, 77 veto 


1.5 to 3.0 


Track finding efficiency 


2.6 to 5.9 


Kaon FID 


n 
O.U 


D° mass 


2.3 to 4.4 


O* -D° mass difference 


2.5 to 6.7 


Dalitz structure 


0.0 to 5.0 


Helicity angle 9* H 


0.0 to 3.8 


Thrust angle 8^ 


5.5 to 7.3 


B° angle 8* B 


3.0 to 3.8 


AE 


8.6 to 12.0 


mES 


2.3 to 4.0 


Simulation statistics 


2.0 to 4.7 


Sum 


16.5 to 19.8 



simulation, since the AS-toes distributions of several 
categories of BB background peak inside the signal box. 
After counting the MC events and scaling the number to 
79.9fb~ 1 , a total of 9.4 ± 1.7 background events is ex- 
pected for all six modes combined. Of those, 2.9 events 
originate from B° — ► D*°ir°, 5.1 events from other BB 
decays, and 1.4 events from qq events. The breakdown 
for each channel is given in Table ITT1 

The estimate of the number of background events is 
cross-checked by two studies, one based on events in the 
grand sideband, and the other based on events in the 
signal box using a control sample of £>*°7r° events. The 
first study results in ratios of data-to-MC events ranging 
from 1.0±0.3 to 1.5±0.2 for the various D*° decay modes, 
and a ratio of 1.2±0.1 for all modes combined. Taking the 
uncertainties into account, data and MC simulation do 
not disagree significantly. For the second study, B° — > 
D*°tt° events are selected by loosening some selection 
requirements and by inverting the 7r veto: we now keep 
events in which a photon combined with the hard photon 
forms a reasonable 7r° candidate. The number of events 
seen in the signal box is usually found to be lower in 
data than in MC simulation with data-to-MC ratios from 
0.3±0.3 to 1.2±0.7 for the various D*° decay modes and 
0.6 ± 0.2 for all modes combined. 

We observe 13 events in the signal box. Figure [31 
presents the AE and 7ties distributions with all selec- 
tion requirements applied. The Monte Carlo simulation 
is shown with separate contributions from B Q — > D*°7r°, 
other BB, and qq events. 

The branching fractions are determined in a 
frequentist-model approach, modified based on Ref. [lflj . 
Besides taking the systematic uncertainty in the effi- 
ciency and the statistical uncertainty in the background 
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TABLE II: Results for individual modes and all modes combined. The upper limit is given for 90% C.L. 



Branching fraction Relative systematic 
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D° -> K--K+ 
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16.5 
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19.8 
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0.0 to 1.3 
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D° -» if-7r + 7r + 7r" 


4.6 


17.3 


2.0 ±0.3 


0.7 ±0.1 


0.5 to 2.0 
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3.9 


D*° -> B° 7 
















D° -> A"-7r+ 


1.4 


17.3 


3.8 ±0.7 


1.6 ±0.4 


0.4 to 1.6 


2 


8.0 


D° -» A-^+tt 


4.9 


19.6 


0.9 ±0.2 


2.4 ± 1.2 


0.1 to 1.2 


3 


14.8 


£>° -» A"-7r + 7r + 7r- 


2.8 


17.7 


1.7 ±0.3 


1.2 ±0.2 


0.2 to 1.7 


5 


20.3 


All modes combined 


23.9 


16.8 


1.8 ±0.3 


9.4 ± 1.7 


0.4 to 1.3 


13 


2.5 



estimate into account, the background expectation value 
is also shifted by a factor selected from a flat distribu- 
tion of the range determined by the data-to-Monte Carlo 
ratios (see Tabic [nj. When combining all six modes, 
this shift comes from the range 0.4 to 1.3 (derived from 
0.6 ±0.2 and 1.2 ±0.1) and is applied coherently for each 
of the modes. Wc assume that 50% of the T(4S) mesons 
decay into neutral BB pairs. Figure 0] displays 1— C.L. 
versus the assumed branching fraction. The significance 
of this measurement, i.e., 1— C.L. at branching fraction 
zero, is 0.86. The central value of the branching fraction 
of B° -> £>*°7 is (1.0iJ]g) x 10" 5 , which is consistent 
with zero. The upper limit on the branching fraction is 
B(B° -> £>*°7) < 2.5 x 10~ 5 at 90% confidence level and 
is in agreement with the theoretical expectations. 




AE(GeV) m ES (GeV/c 2 ) 



FIG. 3: AE (left) and toes (right) distributions for data 
(points) and MC simulation (shaded histograms). All selec- 
tion requirements are applied including the toes signal box 
requirement for the left plot and the AE signal box require- 
ment for the right plot. 
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FIG. 4: 1— confidence level versus the assumed branching frac- 
tion. The shaded areas are the 68% and 95% probability re- 
gions. The 90% C.L. is marked with an arrow. 
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